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Abstract: In this paper, we study the systematics of the 2f states in the N = 82 even-even isotones with proton 
numbers between 52 and 72. We calculate the level energies of the 0*, 27 states and the electric quadrupole reduced 
transition probabilities B(E2; 2; 一 07% in the framework of the nuclear shell model with a monopole- and multipole- 
optimized realistic interaction. Our calculations yield good agreement with the experimental data and show a 2.5 MeV 
gap at Z = 64 subshell closure in '“°Gd. We predict that the B(E2;2} — 07) value for '*°Gd is close to those for '“*Nd 


and '4Sm, and the values increase rapidly from ‘sDy to !7Yb. 
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1 Introduction 


The N = 82 isotones are one of the longest chains of 
neutron semimagic nuclei in nuclide chart. The systemat- 
ics of the 27 states in the N = 82 even-even isotones with 
Z > 50 has attracted much interest. Previous experiments 
measured excitation energies of the 27 states [denoted by 
E,,(27) in this paper] for '**Te-'*Hf, and showed E,(2}) in 
146 Gd is much larger than those in any other N = 82 isotones 
= indicating the proton Z = 64 subshell closure for nuclei 
in this region "Asa comparison, in the proton semimagic 
Sn isotopes with 50 < N < 82, the largest E, (27) is at "Sn, 
and E,(2*) in ''*Sn is the second largest. 

Another important observable for nuclear properties of 
low-lying states is the electric quadrupole reduced transi- 
tion probabilities B(E2;2; 一 07). For example, in the Sn 
isotopes, there is a particular measurement of interest: the 
B(E2;2; 一 07) values, which show a shallow minimum 
at ''6Sn. This minimum has sparked significant attention 
among researchers em and was explained in two different 
ways: a result of a possible soft neutron N = 64 subshell 


[10] ‘ i A 
closure , or due to proton excitation from the Ogon orbit 
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"" For the N = 82 isotones, the measured B(E2; 27 > 07) 
values first increase from "Te to !*°Ce and then stay nearly 
constant up to 144Sm. Unfortunately, no measurements ex- 
ist for isotones with A > 144. This raises the question of 
whether a shallow minimum in B(E2;2; 一 07) exists at 
Z = 64 for the N = 82 isotones. 


In order to comprehensively understand the evolution 
of proton shell structure in N = 82 isotones, theoretical ef- 
forts have been devoted to the systematical study of the 2 了 
states. For example, Holt et al. na performed shell model and 
QRPA calculations. Their results showed that, using a real- 
istic interaction derived from the meson-exchange potential, 
the shell model predicted a slight increase in the E, (27) val- 
ues from !**Te to '4°Gd. On the other hand, the QRPA cal- 
culation produced a constant E,(2}) value. Coraggio et al. 
= performed shell model calculations using a realistic inter- 
action derived from the CD-Bonn potential. They showed a 


maximum E,(2f) value at '**Dy. 


The purpose of this paper is to study the systematics of 
the 27 states forthe N = 82 even-even isotones, in the frame- 
work of the shell model. We calculate the E, (27) values and 
B(E2; 2; — 07) transition strengths using a monopole- and 
multipole-optimized effective interaction that was recently 
developed in our previous work "This paper is organized 
as follows. In Sec. 2 we give a brief introduction to the 
framework of shell model, the Hartree-Fock, and the parti- 
cle number conserved BCS. In Sec. 3 we present our cal- 
culation results including the E (27) values, relative binding 
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energies, and B(E2; 2; — 07) transition strengths. We show i.e., 
our results are in good agreement with the experimental data. Less 
|p) = —(P})” |0), 3) 
In Sec. 4 we summarize our results. XN 


2 Framework 


We perform full shell model (SM) calculations for 34Te, 
136Xe，138Ba，l140Ce，142Nd，I4Sm，146Gd，148Dy，150Er， 
152Yb, and 154Hf with valence protons outside doubly magic 
core '*Sn in the 0g7p1dsp 1d3n25 „Oh, ,n shell using the 
BIGSTICK code’. The maximum dimension of the SM 
calculation reaches 1.6 x 10” at '*?Dy. In our previous work 
u we derived a monopole- and multipole-optimized effec- 
tive interaction based on the realistic JIS6PNA interaction 
I The Coulomb interaction between the protons is included 
in this interaction, similar to previous studies in Refs. [12] 
and [13]. This effective interaction reproduces the binding 
energies, low-lying level energies, the electric quadrupole 
moments, the electric quadrupole reduced transition prob- 
abilities, the magnetic dipole moments, and the magnetic 
dipole reduced transition probabilities for both the even-even 
and odd-mass nuclei with N = 82. In this work, we use this 
interaction to calculate o- and 27 -state level energies and the 
B(E2; 27 
transition operator is defined by T(E2) = e,?”Y,, where e, 


— 07) transition strength. The electric quadrupole 


is the effective charge of valence protons. We take e, = 1.6 
and the empirical formula of the harmonic oscillator length 
b = (1.012A"3)!? fm. 

For comparison, we also calculate the ground state en- 
ergies in the Hartree-Fock (HF) and the number conserved 
BCS (NBcs) ™” 
start with the HF calculation in the Og- 1dsp 1d3n25 0A4 1 


using the same effective interaction. We 


shell with Kramers degeneracy, i.e., our HF calculation pro- 
duces time reversal single-particle partners without enforc- 
ing additional constrains such as shape and orientation. The 
HF single-particle state obtained from the calculation can be 
written as a transformation of the SM single-particle states: 


ee (1) 


Here U is the transformation matrix. 
The building blocks of the NBCS are collective pairs in 
the HF basis, i.e., 


pi = 


MD 1 一 


> oie, (2) 
a 


where v, is the pair structure coefficient. In the NBCS, the 
ground state of 2N valence protons is an N -pair condensate, 


where yy is the normalization factor. It is worth empha- 
sizing that the particle number is exactly conserved in the 
NBCS state described in Eq. (3), which has similarities to 
the seniority-zero state in the generalized seniority scheme 
2l Unlike traditional BCS calculations in nuclear structure 
theory, the NBCS state does not require a numerical particle 
number projection. 

The Hamiltonian of NBCS for valence protons can be 


written as 


H= 2: Eap 人 2 人 5 + 
(4) 


where €,, and V,,,; are single-particle energies and two- 


body matrix elements, respectively, in the HF single-particle 
basis. The energy of the NBCS state described in Eq. (3) can 
be written by 


E=(d| A |) 


[o] 
= F Qena + Vaaa) (: = a ) 
a XN 

oz 天 1 [ap] 5) 
Xn-1 ( 

+ VoaggN Vat 

2 | BB B tn 


[Al _ [ap] 
x + lt xx 
XN 


Here }), represents the summation over a or &, where a 
[ap] 3 


+Vi gap + Vaßaß) (1- 


and & are degenerate time-reversed pairs. we 1 and XN 
æ- and af-orbit blocked normalization factors, oe 
The pair structure coefficient v, is determined by minimiz- 
ing the energy of the NBCS state. The formulae for the vari- 
ational principle of the NBCS can be found in Refs. [18-19] 
and were extended to the case of open-shell nuclei in Ref. 
[20]. Since the N = 82 even-even isotones are nearly spher- 
ical nuclei, the result of the NBCS is very close to that of the 


generalized seniority scheme with seniority zero. 


3 Results 


Table 1 compares for the relative nuclear binding ener- 
gies (the binding energy difference between the N = 82 nu- 
clei and 1? Sn) from the experimental data, our full SM, HF, 
and NBCS. The relative binding energy data are derived by 
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Table 1 Nuclear binding energy difference (in MeV) between the 10 
N = 82 nuclei and 32Sn. The theoretical values are ob- (a) Ogın 
tained by our calculations of full SM, HF, and NBCS with Sp 
the effective interaction in Ref. [14]. 5 L L | | 
20 
Nuclide Expt. SM HF NBCS 5 st (b) 1dsn 
Te 20.548 20.552 19.319 20.548 8 
136Xe 39.004 39.019 36.985 38.993 3 i | | | | | | | 
Ba 55.396 55.421 53.016 55.363 3 o i 
Ce 69.767 69.792 67.410 69.697 © | (e) Oh BOF 
INd 82.200 82.228 78.748 82.079 5 
48m 92.771 92.846 89.434 92.617 
46Gd 101.446 101.568 99.300 101.275 on Ga EI 5 i 154 
8Dy 107.774 107.936 104.482 107.645 A 
SOEr 112.298 112.403 108.553 112.116 
152Yb 115.289 115.116 111.305 114.878 Fig.2 The expectation value of the occupation number of the 
154Hf 116.287 116.158 112.995 115.972 087, 1d;n, and Oh, |» orbits in the OF and 27 states. 


subtracting the electron binding energy from the atomic bind- 
[24] 
. The 


relative binding energies obtained by our SM calculation are 


ing energy value compiled in the AME2020 database 


very close to the data, with a root-mean-square deviation of 
98 keV. The relative binding energies obtained by the NBCS 
are in good agreement with the data and the SM result, but 
those obtained by HF are 1-4 MeV smaller. This indicates 


the pair correlation is very important here. 
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Fig. 1 The excitation energies of the 2f states, E,(27), for the 


‘SM-Holt” 


and “SM-Coraggio” represent the shell-model results in 


N = 82 even-even isotones with Z > 50. 


Refs. [12] and [13], respectively. The experimental data 
are taken from the ENSDF database °” 


Fig. 1 compares for the E,(27) values from the ex- 
perimental data, the SM results in this work and in Refs. 
[12-13]. The data are taken from the ENSDF database = 
We see for the entire range of the N = 82 even-even iso- 


tones, the E, (27) values obtained by our SM calculation are 


in good agreement with the data, with a root-mean-square 
deviation of only 49 keV. The experimental data exhibit two 
distinct peaks of E,Q7), with the highest peak observed at 
'46Gd with an energy of 1.972 MeV, and the second peak at 
140Ce with 1.596 MeV. Our SM calculation successfully re- 
produces this phenomenon, with our values for the two peaks 
being 1.961 MeV and 1.664 MeV, respectively. On the other 
hand, the SM calculation in Ref. [13] predicted E (27) in 
good agreement with the data for '**Te-'**Sm, but it failed 
to reproduce the maximum peak at '°Gd: the theoretical 
E,(27) value for !4°Gd is 0.3 MeV smaller than the data, and 
those for !*?Dy, °Er, and !?Yb are ~ 0.2 MeV larger than 
the data. The effective interaction used in Ref. [13] does 
not well reproduce the evolution of the proton shell structure. 
The calculation in Ref. [12] predicted E,(27) values increase 
slightly from 134Teto '4°Gd, but the results for A > 146 were 


not reported. 


The expectation value of the occupation number of the 
Ogn, 1dsp, 1d3p, 2817, and Oh, 1,7 orbits is obtained by our 
SM calculation. The occupation numbers of the 1d3,. and 
25,/. orbits are very small in the Of and 2} states of !**Te- 
4h f In Fig. 2 we see that the occupation numbers of the 
0g7 and 1ds,, orbits increase rapidly and approach the satu- 
ration values as A increases, while the number of the 0h, ın 
orbit almost vanishes for A < 146 and subsequently in- 
creases. The occupation number of the 1ds orbit in the 2 
state of '4°Gd is 0.7 smaller than that in the ground state, 
whereas the occupation number of the Oh, |„ orbit in the 2 了 
state is 0.7 larger than that in the ground state. The 2 state 
n \46Gd is conventionally considered a one-phonon excita- 
tion of the 07 ground state or the generalized seniority two 
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"" The motion of the valence proton from the 1d,,, to the 


Oh,» orbit plays a role in the formation of the 27 state. 

In order to study the proton shell structure for the N = 
82 isotones, we derive the HF single-particle energy for the 
SM orbit. The procedure is as follows. Using Eq. (1) we 
have 


Ej m, = > (Oya) Ea» (6) 


where e, is the calculated single-particle energy of the HF 
orbit, and sj m, is the single-particle energy of the SM orbit. 
Our calculation shows that for the nearly spherical even-even 
nucleus with N = 82, the € Jm values with the same j, but 
different m, are close to each other. Based on this obser- 
vation, we calculate the average sj „, value for each j, but 
differing m,: 

E 三 : BIER (7) 

a 2j, + 1 JaMa 


Ma 


We assume that €;, represents the HF single-particle energy 
of the j, orbit. 


Eja (MeV) 


1 1 1 1 1 1 L 1 1 1 1 
134 138 142 146 150 154 
A 


Fig.3 The evaluated HF single-particle energy €, of the 0g7,, 


lds, 1d3)9, 2819, and Oh11 orbits. 


Fig. 3 shows the results of the HF single-particle ener- 
gies of the 0g7, 1dsp, 1d3p, 2512, and Oh} |, orbits in the 
N = 82 even-even isotones. We see a pronounced increase 
in the energies across the isotones 134Te-154Hf with values 
ranging from —9 to —6 MeV at lower energies and increas- 
ing to —3 to 2 MeV at higher energies. This increase arises 
due to the repulsive monopole interaction between identical 
valence nucleons. We find a large gap between the 1d,,, and 
OA 1 orbits in Gd: the difference between £; dyp and Eon 
is equal to 2.5 MeV. This large gap produces the Z = 64 sub- 


shell effect, which manifests as a very large excitation energy 
ofthe 27 state. Interestingly, this gap gradually decreases as 
nuclei move away from 146Gd, revealing a trend in the evo- 
lution of the proton shell structure. Similarly, we find a 1.6 
MeV gap between the lgyp and 1ds,, orbits in !*°Ce, which 
is responsible for the large excitation energy of the 27 state. 


0.08 - 
se 0.06- 
> 
5 a 
A 0.04 - “7 von Garrel06 
— = SM-This work $ Chenl7 
2 r =-= SM-Holt $ Eccleshall66 
0.02- =." SM-Coraggio f Kim92 
á $ Danchev1l $ Pitthan73 
$ Stuchbery13 $ Kindleben78 
0.00 1 1 1 1 1 1 1 1 1 1 1 
134 138 142 146 150 154 
A 
Fig.4 The electric quadrupole reduced transition probabilities 


B(E2;2} — 07) for the N = 82 even-even isotones with 
Z > 50. ‘SM-Holt” and “SM-Coraggio” represent the SM 
results in Refs. [12] and [13], respectively. The experi- 
mental data are taken from Refs. [26-34]. 


Fig. 4 compares for the B(E2; 2f 一 07) values from 
the data, the SM results in this work and in Refs. [12-13]. 
The data are taken from Refs. [26-34]. We see for the lighter 
N = 82 nuclei with A < 144, our SM result is in good 
agreement with the data within the experimental uncertain- 
ties, except that for 138Ba the B(E2) value obtained by our 
calculation is 25% smaller than the data. On the other hand, 
the SM results in Refs. [12-13] for "Te, '?°Xe, '?®Ba, and 
140 Ce are in good agreement with the data, but those for !*?Nd 
and '**Sm are too large. The B(E2;2} — Of) value for 
the isotones heavier than ‘Sm has not been measured. Our 
SM calculation predicts that the B(E2;27 — Of) value for 
'46 Gd is close to those for '47Nd and '**Sm, and the value in- 
creases rapidly from !*®Dy to !* Yb and drops at !*Hf. The 
Z = 64subshell closure does not lead to a shallow minimum 
of B(E2;2} > Of) at !*%Gd. This feature is different from 
the feature in Sn isotopes that there is a shallow minimum 
near N = 64. 
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w 
= 


Q. Y. Chen et al: Systematics of 27 states in N = 82 even-even isotones ESS 


4 Summary 


In this paper we calculate low-lying level energies and 
B(E2; 27 一 07) values for N = 82 even-even isotones in 
the shell model with an effective interaction. The experi- 
mental data exhibit two distinct peaks of E, (27) at Gd and 
140Ce, Our shell model results for E,(27) and binding ener- 
gies are in good agreement with the experimental data, while 
previous shell model results are not as good as ours. Our cal- 
culation also well reproduces the B(E2; 27 — 07) for 134Te, 
136Xe, 140Ce, Nd, and ‘Sm. The calculation and the ef- 
fective interaction used in this work describe correctly the 
proton shell structure for the N = 82 isotones. 


Our HF calculation shows a 2.5 MeV energy gap be- 
tween the Ids) and Oh, ,„ orbits in '*°Gd. This large gap 
results in the formation of the Z = 64 subshell closure. Sim- 
ilarly, we find a 1.6 MeV energy gap between the Og- and 
Ids) orbits in '4°Ce. Our shell model calculation predicts 
the B(E2; 2, 一 0f) value for '“°Gd is close to those for 
142Nd and !“ Sm, and the values increase rapidly from !*?Dy 
to "Yb. The existence of the Z = 64 subshell does not lead 
to a shallow minimum of B(E2; 2} > 07) at “Gd. 
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N = 82 偶偶 同 中 子 素 27 态 的 系统 性 


陈 青云 mi, NEAR, Mae, EPA 
(1. 同济 大 学 物理 科学 与 工程 学 院 , 上 海 200092; ) 


BE: 在 本 文中 我 们 研究 质子 数 在 52 和 72 之 间 的 N = 82 偶偶 同 中 子 素 2+ 态 的 系统 性 。 我 们 用 原子 核 壳 模 
型 结合 单 极 和 多 极 修正 的 真实 相互 作用 计算 了 OF 和 21+ 态 能 级 以 及 电 四 极 约 化 跃迁 几率 BE 27 > 07) R 
们 的 计算 结果 与 实验 数据 非常 吻合 。 我 们 的 研究 表明 在 “Gd 中 Z = 64 子 壳 的 能 辽 为 2.5 MeV。 我 们 预测 
“6Gd 的 B(E2;2* 一 01) 值 与 INd 和 JSm 的 值 接近 ， 从 1%Dy 到 52Yb 的 B(E2) 值 迅速 增 大 。 

关键 词 : WRAN = 82 同 中 子 素 ;21 态 能 量 ; 电 四 极 约 化 跃迁 几率 
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